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SUMMARY 


The theory of anisotropic elasticity was used to evaluate the anisotropic stress concentration factors of a composite 
laminated plate containing a small circular hole. This advanced composite material was used to manufacture the 
X-29A forward-swept wing. Observe that the usual isotropic material stress concentration factor is three. However, 
for composite material, it was found that the anisotropic stress concentration factor is no longer a constant, and 
that the locations of maximum tangential stress points could shift by changing the fiber orientation with respect to 
the loading axis. The analysis showed that through the lamination process, the stress concentration factor could be 
reduced drastically, and therefore the structural performance could be improved. Both the mixture rule approach 
and the constant strain approach were used to calculate the stress concentration factor. The results predicted by the 
mixture rule approach were about twenty percent deviate from the experimental data. However, the results predicted 
by the constant strain approach matched the testing data very well. This showed the importance of the inplane shear 
effect on the evaluation of stress concentration factor for the X-29A composite plate. 
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INTRODUCTION 

The X-29A advanced technology demonstrator is sponsored by the Defense Advanced Research Projects Agency 
with support from NASA and the Air Force. 

The X-29A features a unique forward-swept wing (Fig. 1 ), made of composite materials which offers weight 
reduction of as much as 20 percent in comparison with convention aft-swept wings. 

A forward-swept wing is prone to structural divergence, because as dynamic pressure increases, forces tend 
to bend the leading edge up. If a divergent speed were reached, a cycle of leading edge bending, increased local 
angle of attack and greater wing loading could grow to cause a structural failure. The wing’s high rigidity prevents 
divergence from occurring within the X-29A’s flight envelope. 

Because of their importance in aircraft design application, laminated, continuous-fiber reinforced-resin matrix 
composites containing through cutouts have been the subject of considerable study (refs. 1 to 5) in this paper, 
anisotropic plate theory was used to calculate the anisotropic stess concentration factors (SCF) for the X-29A com- 
posite plate containing a circular hole. 

NOMENCLATURE 

Aij equivalent modulus or extensional stiffness for a multidirectional laminate 

E\ modulus of elasticity of anisotropic plate in axis-1 direction 

E 2 modulus of elasticity of anisotropic plate in axis-2 direction 

Ei modulus of elasticity of laminated composite plate in axis-1 direction 
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V\2,V2\ 

5)2,521 

VLT 

Ci 
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modulus of elasticity of laminated composite plate in axis-2 direction 

transformed modulus of elasticity of single-ply in axis-1 direction 

transformed modulus of elasticity of single-ply in axis-2 direction 

modulus of elasticity of anisotropic plate in a direction 

modulus of elasticity of single-ply to fiber direction 

modulus of elasticity of single-ply transverse in fiber direction 

shear modulus of anisotropic plate associated with {1,2} coordinate system 

shear modulus of laminated composite plate associated with {1,2} coordinate system 

transformed shear modulus of single-ply in {1 , 2 } coordinate system 

shear modulus of single-ply associated with {L, T} coordinate system 

stress concentration factor 



stress resultant 

reduced stiffness of single-ply 

transformed reduced stiffness of single-ply associated with {1,2} coordinate system 
Poission’s ratio of anisotropic plate associated with {1,2} coordinate system 
Poission’s ratios of laminated composite plate associated with {1,2} coordinate system 
transformed Poission’s ratios of single-ply in {1 , 2 } coordinate system 
Poission’s ratio of single-ply associated with {L, T} coordinate system 
constant strain 
stress in a direction 
remote tensile stress 


ANALYSIS 

Let axes 1,2 be the principal coordinate axes of the laminated plate, and let axes L,T be the principal material axes 
of the single composite ply shown in figure 2. 

For an anisotropic plate containing a circular hole subjected to remote uniaxial tensile stress acting at an 
angle <j) with respect to the principal elastic axis 1 of the plate (fig. 3), the tangential stress, o a and tangential stress 
concentration factor, K = o a /a^ along the circular hole boundary may be expressed as (ref. 3) 

_ n 

K = — = -=r{[ — cos 2 <j) + (k + n) sin 2 (j)]k cos 2 a+ [( 1 + n) cos 2 0 — A: sin 2 <f>] sin 2 a (1) 

- n( 1 + k + n) sin <}> cos (j> sin a cos a} 

where E a is the modulus of elasticity in the a direction (fig. 3) given by 
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where k and n are defined by 


k = ~H\H 2 
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where i = \/— 1 , and and H2 are the complex roots of the anisotropic plate characteristic equation 
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For isotropic materials k = 1 and n= 2 , and the stress concentration factor K (eq. ( 1 )) reduces to 


K = — =1—2 cos 2 (a — <f>) (6) 

^OO 

which gives K = - 1 at ( a - <f>) = 0 or 7r, and if = 3 at ( a — </>) = ±7r/2 . 

To evaluate the modulus of elasticity of a laminated plate, both the mixture rule approach and the constant strain 
approach could be used. In the mixture rule approach, the transformed ply-elastic constants {Ei , Ez , G\z , ^12 , vi \ } 
with respect to the {1 , 2} system can be related to the material constants {El,Et,Glt,vlt, vtl} with respect to 
the {L, T} system through the following equations (refs. 6 and 7 ). 
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If the composite plate is made of iV number of single plies with different fiber orientations, then by using the mixture 
rule, the engineering elastic constants {ij| , Ez , G\z , v\2 , ^21 } for the composite plate can be written as 



( 8 ) 
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In the constant strain approach, it is assumed that the strain remains constant across the laminate thickness and 
the inplane stress-strain relation for a laminate is used and it is actually the stress resultant versus inplane strain 
relation. 


N\ = All ei + Ai2 C2 + A\6 £6 

N2 - A 21 C| + Anti + ^26 C6 (9) 

Ne = A$l £| + A62€l + A(& €6 

where Aij are defined by (refs. 6 and 7) 

N 

Aij = J2(QiMzk-z k - l ) *,; = 1 , 2,6 ( 10 ) 

*=1 

in which 


Qu = Qu cos 4 0 + 2(Q, 2 + 2Q66) sin 2 0cos 2 6 + Q 22 sin 4 9 
Q \2 = (Qu + Q 22 — 4Q66) sin 2 9 cos 2 0 + Qi 2 (sin 4 9+ cos 4 9) 

Q 22 = Qu sin 4 0+ 2( Q 12 + 2Q 66 ) sin 2 0cos 2 0+ Q 22 cos 4 0 
0i6 = (Qu - Q 12 - 2Q66) sin 0COS 3 0+ (Q 12 - Q 22 + 2Qa;) sin 3 0cos 0 
Q 26 = (Qu - Q 12 - 2Q66) sin 3 0cos 0+ (Q I2 - Q 22 + IQes) sin 0cos 3 0 
Q66 = (Qu + Q 22 - 2 Qi 2 - 2Q 66 ) sin 2 0cos 2 0+ Q66(sin 4 0 + cos 4 0) 
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( 11 ) 


( 12 ) 


The calculation of the effective engineering elastic constants {jji , Ei , G\i , £12 } is performed by relating the 
compliance components to inplane engineering constants under uniaxial tension along the 1-axis (ref. 8). 
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RESULTS 


The X-29A forward-swept wing composite plate is made up of 40-plies with the total thickness of .56 cm (.22 in). 
The stacking sequence and the ply-engineering elastic constants are given by 

[ ±45 10 4 1 ± 45 190 2 1 ± 45 10 4 1 ± 45 |0 2 ] s 

El = 18.76 x 10 6 psi 
Et = 1 .57 x 10 6 psi 
Glt = 0.82 x 10 6 pst 
vlt = 0 .312 


By replacing {E\ , E 2 , Gn , un }, respectively, with {Ei , E 2 , G 12 , v\i } in equations (1), (2), (3), and (4), the 
tangential stresses a Q around a circular hole in a laminated X-29A composite plate were calculated for three loading 
cases: <f> = 0 (loading in axis-1 direction), <f>= f, and 0 = j (loading in axis-2 direction). The results obtained from 
the constant strain approach and the mixture rule approach are ploted in figures 4 to 6 and figures 7 to 9 respectively. 
Figure 4 shows the plot of a Q for the laminated X-29A composite material when the plate is under uniaxial tension 
in the composite elastic axis-1 direction. The maximum stress concentration factor K for the laminated X-29A 
composite plate reached the peak value of 3.614 (larger than 3) at two locations (a = 90° and a = 270°). When the 
loading axis is <f> = f oblique to the composite axis-1 (fig. 5), the maximum stress concentration factor K reaches the 
value of 3.033 at two locations (a= 1 20 0 and a = 300 °). When the loading axis is parallel to the composite elastic 
axis-2 (fig. 6), the maximum stress concentration factor K drops to the value of 2.708 at two locations (a = 0° and 
a = 180°). 

Again, figure 7 shows the plot of a a for the laminated X-29A composite material when the plate is under uniaxial 
tension in the composite elastic axis-1 direction. The maximum stress concentration factor K for the laminated X- 
29A composite plate reached the peak value of 4.475 (higher than 3.614 in figure 4) at two locations (a = 90° 
and a = 270°). When the loading axis is <f> = f oblique to the composite axis-1 (fig. 8), the maximum stress 
concentration factor K drops to the value of 2.977 (lower than 3.033 in figure 5) at two locations (a = 100° and 
a = 290 °). When the loading axis is parallel to the composite elastic axis-2 (fig. 9), the maximum stess concentration 
factor K reaches the value of 3.021 (higher than 2.708 in figure 6) at two locations (a = 0° and 180°). Finally, 
figures 4 to 6 and figures 7 to 9 are summarized in Table 1 and Table 2, respectively. 

For comparison purposes, similar calculations were made for a single-ply of the X-29A composite using equation 
(1). When loading is along the fiber direction (axis L, fig. 10), the maximum stress concentration factor K reaches 
the peak value of 6.401 at two locations (a = 90° and a = 270°). When the loading is f oblique to the fiber direction 
(fig. 11), the maximum K value is 3.988 at two locations (a = 100° and a = 280°). When the loading direction 
is transverse to the fiber direction (axis T, fig. 12), the maximum value of K is -3.457 at two points (a = 90° 
and a = 270 °). The negative sign of K value represents a compressive stress concentration. Table 3 is used to 
summarize the results from figures 10 to 12. For reference purposes, the stress concentration factor of isotropic 
material is plotted in figures 13 and 14. 

The stress concentration factors evaluated from different approaches (mixture rule and constant strain) discussed 
previously were compared by performing simple coupon tests. The width W of the rectangular specimen is 3.81 
cm (1.5 in) and the diameter of the small central circular hole is .635 cm (0.25 in). The comparison of stress 
concentration factors between theoretical predictions and experimental results are listed in table 4. The comparison 
of stress concentration factors between single-ply and laminated plate are shown in table 5. 

A simple FORTRAN program for calculating anisotropic stress concentration factors is listed in the appendix. 
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CONCLUSION 


The theory of anisotropic elasticity was used to evaluate the anisotropic stress concentration factors for single-ply 
and laminated OBX-29A (forward-swept wing) research aircraft composite plates, each of which contained a small 
circular hole. 

It is well known that the usual isotropic material stress concentration factor is three. However, the analysis 
showed that the anisotropic stress concentration factor could be greater or less than three for composite materials, 
and the locations of the maximum tangential stress points could shift by the change of fiber orientation with respect 
to the loading axis. 

It was found that through the lamination process the stress concentration factor could be reduced drastically, and 
therefore the structural performance could be improved. The next logical step in the study of anisotropic stress con- 
centration problem would be to know the optimum lamination process to obtaining the minimum stress concentration 
factors of laminate plates. This is a subject that may need further studies. 

Both the mixture rule approach and the constant strain approach were used to calculate stress concentration 
factors. The results obtained by the mixture rule approach were about twenty percent deviate from the experimental 
data. However, the results predicted by the constant strain approach matched the testing data very well. This showed 
the importance of the inplane shear effect on the evaluation of stress concentration factors for the laminated X-29A 
composite plate. A further investigation about the inplane shear effect will need a three dimensional model from 
anisotropic elasticity plus the interlaminar stress analysis. 

The anisotropic stress concentration of laminated plates is a difficult and complicated problem. To obtain a better 
understanding of this physical phenomenon, consideration of the hole size effect and utilization, the theory of linear 
elastic fracture mechanics, and the theory of micromechanics is imperative. 
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TABLE 1 . Stress Concentration Factors 
for an X-29A Laminated Plate 40-Plies 
[ ±45 10 4 1 ± 45 190 2 1 ± 45 10 4 1 ± 45 10 2 ], 
Constant Strain Approach 



0° 

45° 

90° 

K 

3.614 

3.033 

2.708 

a 

90°, 270° 

120°, 300° 0°, 180° 


<f>: orientation of loading axis 
K : stress concentration factor 
a: locations of peak stress 


TABLE 2. Stress Concentration Factors 
for an X-29A Laminated Plate 40-Plies 


[ ±45 10 4 1 ± 45 190 2 1 ± 45 10 4 1 ± 45 10 2 ] 4 
Mixture Rule Approach 


<t> 

0° 

45° 

90° 

K 

4.475 

2.977 

3.021 

a 

90°, 270° 

110°, 290° 

0°, 180° 


<f>\ orientation of loading axis 
K : stress concentration factor 
a: locations of peak stress 


Table 3. Stress Concentration Factors for a 
Single-Ply of X-29A Composite 



0° 

45° 

O 

O 

On 

K 

6.401 

3.988 

-3.457 

a 

90°, 270° 

100°, 280° 

90°, 270° 


4>: orientation of loading axis 
K: stress concentration factor 
a: locations of peak stress 
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TABLE 4. Comparison of Stress Concentration Factors 


Loading 

(lb) 

Constant 

strain 

Mixture 

rule 

Experimental 

results 


3.614 

4.475 

3.627 


3.614 

4.475 

3.506 


3.614 

4.475 

3.567 


3.614 

4.475 

3.546 

WSM 

3.614 

4.475 

3.506 


<j> = 0°: orientation of loading axis 
a = 90°, 270°: locations of peak stress 


TABLE 5. Comparison of Stress 
Concentration Factors 


Single 

Laminated plate 

ply 

Constant strain Mixture rule 

6.401 

3.614 4.475 


0 = 0°: orientation of loading axis 
a = 90 °, 270°: locations of peak stress 
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APPENDIX 

C NSCF.FOR NASA X29 WING STRESS CONCENTRATION FACTOR 
WRITE (*,10) 

10 FORMAT ('INPUT El, E2, V12, G12 f , / ) 

READ (*,20) El , E2, V12, 612 

20 FORMAT ( 4E12.5) 

WRITE (*,21) 

21 FORMAT (2X , ' El ' , 4X , » E2' , 4X , * V12' , 4X , » G12' ) 
WRITE (*, 22) El , E2, V12, 812 

22 FORMAT (4E12.5) 

E0E=E1/E2 
E0G=E1/G12 
SK=CEOE) **. 5 

SN= ( 2* ( SK-V 1 2 ) +EOG ) ** . 5 
PI =3. 14159265359 
PH 1=0. 

DO 100 L=1 , 3 
RPH I =F‘H I *P I / 1 80 . 

SS I =S INC RPH I ) 

CCO=COS C RPH I ) 

IF ( ABS (SSI ) . LE. .0000001) GO TO 101 
IF ( ABS (CCO) . LE . .0000001) GO TO 102 
GO TO 110 

101 SSI=0. 

GO TO 110 

102 CC0=0. 

GO TO 110 

no cc 2 =cco*cco 

SS2=SSI*SSI 

CCSS=CCO*SSI 

AN=0. 

DO 30 K=1 , 73 
ALF A= AN*P I / 1 80 . 

SI=SIN( ALFA) 

CO=COS C ALFA ) 

IF C ABS (SI). LE. .0000001) GO TO 31 
IF ( ABS (CO). LE. .0000001) GO TO 32 
GO TO 40 

31 SI=0. 

GO TO 40 

32 C0=0. 

GO TO 40 
40 CS=SI*CO 
C2=C0*C0 
S2=S I *S I 
C4=C2*C2 
S4=S2*S2 

EAF A=E 1 / C S4+E0E*C4+ ( EOG-2* V 1 2 ) *S2*C2 ) 

AA= (-CC2+ CSK+SN) *SS2) *SK*C2 
BB= C C 1 +SN) *CC2-SK*SS2) *S2 
DD=-SN* ( 1 +SK+SN ) *CS*CCSS 
SCF= (EAFA/E1 ) * C AA+BB+DD) 

WRITE (*,200) 

WRITE (*,210) PHI, AN, SCF 
AN=AN+5. 

WRITE (*,300) 
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300 FORMAT (2X, ' ANGLE ALPHA =AN =LOCATION OF STRESS POINT') 
WRITE (*, 310) AN 
310 FORMAT (E12.4) 

30 CONTINUE 

WRITE (*,350) 

350 FORMAT (2X, 'ANGLE PHI= LOADING AXIS') 

WRITE (*, 360) PHI 
360 FORMAT (2X,E12.4) 

PHI=PHI+45. 

100 CONTINUE 

WRITE (*,200) 

200 FORMAT (2X, 'PHI' 8X,'AN', 8X, 'SCF') 

WRITE (*,210) PHI, AN, SCF 
210 FORMAT (3E14.4) 

STOP 

END 
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C XAIJ.FOR NASA X29 WING IN-PLANE STIFFNESS 
WRITE C*,10) 

10 FORMAT C» INPUT INT-L, El, E2, V12, G12’,/) 

READ C*,20) L, El, E2, 012, G12 
DIMENSION ZC24), TK1 C23) 

20 FORMAT (2X, I4,4E12.5) 

DATA C TK1 Cl), 1 = 1,23) / 45. , -45. , 0. , 45. , -45. , 90. , 45. , -45. , 0. , 
A 45. , —45. , 0. , —45. , 45. , 0. , —45. , 45. , 90. , —45. , 45. , 0. , —45. , 45. / 
DATA CZCI), 1=1,24) / -. 1 1 ,-. 1045, -. 099, -. 077, -. 0715, -. 066, 

A -.055, -.0495, -.044, -.022, -.0185, -.011, .011, .0165, .022, .044, 

A .0495, .055, .066, .0715, .077, .099, . 1045, .11/ 

PI =3. 14159265359 
V21=E2*V12/E1 
DUM=1.-V12*V21 
RU=E1/DUM 
R12=V12*E2/DUM 
R22=E2/DUM 
R66=G12 
A1 1 =0. 

A12=0. 

A 16=0. 

A22=0 . 

A26=0. 

A66=0. 

DO 30 K= 1 , L 
TK=PI*TK 100/180. 

SI=SINCTK) 

CO=COS CTK) 

IF CABSCSI) .LE. .0000001) GO TO 9 
IF CABS (GO) .LE. .0000001) GO TO 19 
GO TO 29 

9 SI=0. 

GO TO 29 

19 C0=0. 

29 CD=SI*CO 

C2=C0*C0 
D2=SI*SI 
C4=C2*C2 
D4=D2*D2 

Q 1 1 =R 1 1 *C4+2 . * C R 1 2+2 . *R66 ) *C2*D2+R22*D4 
Q 1 2= C R 1 1 +R22-4 . *R66 ) *D2*C2+R 1 2# C D4+C4 ) 

Q22=R 1 1 *D4+2 . * C R 1 2+2 . *R66 ) *D2*C2+R22*C4 
Q 1 6= C R 1 1 -R 1 2-2 . #R66 ) *C2*CD+ C R 1 2-R22+2 . *R66 ) *D2*CD 
□26= CR 1 1 -R 1 2-2 . *R66 ) *CD*D2+ C R 1 2-R22+2 . *R66 ) *CD*C2 
066= C R 1 1 +R22-2 . #R 1 2-2 . *R66 ) *C2*D2+R66* C D4+C4 ) 

ZD=Z CK+1 ) — Z CK) 

A11=A11+Q11*ZD 
A12=A12+Q12*ZD 
A22=A22+Q22*ZD 
A16=A16+Q16*ZD 
A26=A26+Q26*ZD 
A66=A66+Q66*ZD 
30 CONTINUE 

WRITE C*, 100) 

100 FORMAT C4X, 'All' ,4X, 'A12' ,4X, »A16» ,4X, » A22' ,4X, 'A26' , 
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A 4X,'A66 f > 

WRITE C*:,110) All, A12, A1&, A22, A2&, A66 
110 FORMAT (6E12.4) 

STOP 

END 


(N 



Figure 2. Rotation from material axis L to laminated Figure 3. Tension at an angle to a principal elastic 

plate axis 1. axis 1 of an anisotropic plate with a circular hole. 
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Figured SC F of X-29 laminated plate 
( constant strain). 


Figure 5. SCF of X-29 laminated plate ( con 
slant strain). 



Figure 6. SCF of X-29 laminated plate ( constant 



Figure 7 . SCF of X-29 laminated plate Figure 8. SCF of X-29 laminated plate (mix- 

(mixture rule). ture rule). ;! 
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Figure 9. SCF of X-29 laminated plate ( mixture Figure 10. SCF ofX-29 single ply 

rule). (cj) = 0). 
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Figure 13. SCF of an isotropic plate. 




Figure 14. SCF of an isotropic plate. 
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